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Abstract 

We make use of a unitarized coupled-channel model to analyze the mass distribution data of 
final states in production processes of X (4260). By analyzing the analytical structures of the decay 
amplitudes, we find that the line shape of Z c (3900) signal is related to the combined effect of a pair 
of near-threshold “shadow” poles and the (DD*)^ thresholds, in which the third-sheet pole might 
provides a dominant contribution. As all the coupled channels effects are tuning off, the trajectories 
of these two poles suggest that the Z c (3900) might originate from the attractive interaction of 
(. D*D *) ± through a long-distance interaction, e.g. 7r-exchange interaction, as a “deuteron-like” 
state. There is no nearby pole structure corresponding to the Z c ( 4025) signal in the (D*/)*) 31 
mass distribution. 
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Dynamical origins of near-threshold structures are under heavy debates since more and 


more such signals are observed in experimental explorations jl]. Among others, the ones 
in the charged channels with heavy quarkonium final states are of specific interest in both 
theoretical and experimental investigations because of their possible exotic nature. If these 
structures, dubbed Z c ’s or Z b ’s, are really produced by resonant states, they contain at least 
four quark components in exotic forms of matter such as hadronic molecules, quark-gluon 
hybrids, tetraquarks, and some others Q. 

Intense experimental searches for near-threshold charged exotic candidates started from 
the observations of Z b ( 10610) and Z b ( 10650) structures by Belle in e + e _ —?• T(5S') —* 
T(nS)7T + 7T~ and T(55) —$■ h b (nP)7T + 7T~ processes by analyzing the final Y^nS)^ or 
hbinP)^ mass distributions jd]. The Z c (3900) was reported by BES in e + e“ —> .A(4260) —> 
jLj], and then Z c (3885) in e + e _ —> (DD*) T ^ j^], Z c (4025) in e + e _ — > ( D*D*) t tt ± [s]. 
Since these signals appear just near some open-flavor thresholds, interpreting them as reso¬ 


nances are doubted by some groups and they suggested the signals be caused by the effects 
of threshold cusps (see references { 7 , dj] for example). In fact, the experimental data could 
be reproduced through rescattering models without introducing new states [ 9 , lo|. How¬ 
ever, this mechanism of pure cusp effects is challenged by Ref. [ll], in which the authors 
criticized the calculations mentioned above for not including the higher-order contributions 
which may lead to significant deviations from the lowest-order result. They presented that 
a numerical calculation incorporating higher-order perturbation series in an elastic channel 
could dynamically generate a bound state right below the threshold. Recently, the author 
of Ref. [8] improved his cusp model and claimed that it is not necessary to introduce poles 
to explain the signals of Z c (3900) and Z c (4025) even after the higher-order contributions 
are included [12]. This debate highlights the importance of determining the properties of a 
near-threshold peak signal and also urges us to find a general method to distinguish the cusp 
effect and the near-threshold-pole effects as a reference for future experimental exploration. 

In this paper, we adopt a coupled-channel model in which the unitarity of the amplitudes 
is respected and the analytical structures of the amplitudes could be easily analyzed. The 
contributions from all the relevant two-body hadron loops can be taken into account and 
all the coupled processes could be studied simultaneously. We apply this scheme to study 
the mass distribution data of different final states, J (D/)*) 1 * 1 , and (D*D*) ± in the 
production processes of e + e _ collisions around 4.26GeV at the same time, and the experi- 
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FIG. 1. A pictorial representation of the iterative equation, where the “i”, “j” , and “fc” denote 
the coupled channels. The cross point represents the bare contact interaction and the solid dots 
represent the full four-point interactions. 


mental data could be reproduced in a perfect quality with fewer parameters than those in 
the experimental line-shape analyses. The numerical results of the best fit demonstrate that 
the Z c (3900) signal is found to be generated by a combined effect of two near-threshold poles 
and the ( y DD*) ± threshold, however, there is no pole corresponding to the Z c (4025) signal 
in the ( y D*D*) ± mass distribution data. This scheme could be easily generalized in other 
experimental analyses to determine the possible properties of near-threshold structures. 

We first explain the basic concepts and the setup of the model briefly. The related exper¬ 
imental information mainly come from production processes, but, due to crossing symmetry, 
the Lorentz invariant amplitude in a decay process A i(A(P) —>■ B{jpi) + C(p 2 ) + D(p 3 )) is the 
analytical continuation of that in a scattering process M.(A(P) + D{k )) —> B(pi ) + C(p 2 ), 
where D represents the antiparticlc of D and the four-momentum k = — p 3 , and vice versa. 
In other words, they are the same function of the Mandelstam variables, s — (pi + P 2 ) 2 , 
t = ( pi + p 3 ) 2 , and u = ( p 2 +P 3 ) 2 , defined in different physical regions. 

We can define a simple factorization form for the two-body scattering partial-wave T 
amplitude in “Argand Units” in a matrix form as 


T = G f EG, 


( 1 ) 


where G = diag{G 1 (s), G 2 (s), ■■■ U_and G n (s) = y/ p n (s)f n (s)9(s — Sth,n ) similar to the 
Unitarized Quark Model (UQM) 13]]. p n (s) is the kinematic factors and is represented 


as p n (s) = 2 k n (s)/y/s where k n (s ) is the three-momenta of two particles in the center-of- 
mass (CM) frame of the n-th channel. f n {s) is a simple gaussian form factor defined as 
fn(s) = exp[— /c(s) 2 /2/cq] where /3 = 1/ko could be regarded as the interaction range of 
two particles in the channel and this form factor also serves to regularize the ultra-violet 
divergence in the latter calculation. 9(s — s t h , n ) is a unit-step function and s t h, n is the CM 
energy square of the n-th threshold. The matrix element represents the full partial-wave 
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amplitude excluding the incoming and outgoing kinematic factors, form factors, and the 
unit-step functions in i —>■ j process. Unlike in the UQM, we will not introduce s-channel 
bare “seed” states in the £ function here which means that the poles of the amplitude in 
our later analysis are all dynamically generated from the coupled-channel interactions. If 
the s-channel contribution is dominant, as in the processes discussed in this paper, which 
involve transitions of heavy quark pairs but no annihilations of them, the £ satisfies an 
iterative equation 


£ = A + AII£, (2) 

where the matrix element A y represents the effective contact interaction strength between 
the 2 -tli and j'-th channels including the bare interaction coefficient, polarization summa¬ 
tion factors, and Clebsh-Gordan coefficients in the couplings of angular momenta. II = 
diag{Ui(s), n 2 (s), • • •} where IU(s) is the hadron loop integral function with intermediate 
states of the 2 -th channel. Eq. (121) could be represented pictorially as in FigO which means 
that contributions from all kinds of bubble chains have been summed up, and it could be 


viewed as a simplification of the Lippmann-Schwinger equation as in Ref. 
The coupled-channel unitarity relation 

(T - T ] )/2i = TT ] 


QC|. 


leads to 


which requires 


2 i 


: ((PEG - G t E t G) = G^GG^G, 


( 3 ) 


( 4 ) 


( 5 ) 


Im£ _1 = -GG ] . 

From Eq.((2]), the inverse matrix of £ is represented as 

£ _1 = A _1 (J - All), (6) 

where / is the unit matrix. Thus, 

Imn = GG\ (7) 

which means that the imaginary part of II n (s) function is ImII n (s) = p n (s)f%(s)9(s — s t h,n )• 
The real part of II n (s) function is determined through a dispersion relation 


Ren n (s) = — V 

7r 


ImII n (s) 


dz, 


( 8 ) 
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where V f denotes the Cauchy principal integration. 

Now, we have a unitarized partial-wave coupled-channel scattering amplitude with all 
higher-order loops summed up, in which no s-channel bare state is introduced since X tJ ’s 
only represent the contact interactions of the two channels. But, the resonances or bound 
states could be produced dynamically, if the interaction strengths are strong enough. The 
poles of the amplitude, evaluated by the zero points of determinant of (/ — All), can appear 
on the real axis of the first Riemann sheet for bound states or in the complex s-plane on the 
other non-physical sheets for resonances. 

The pole structures of the amplitude could be analyzed by analytically continuing the 


amplitude into the complex s-plane as widely discussed in the literature jl3|, [l6|, [l7j • Every 


physical cut will double the number of Riemann sheets in the analytical continuation, which 
means that there are 2 n Riemann sheets for an n-channel amplitude. To End the poles 
on the closest Riemann sheet attached to the physical region, one need to analyze the 
analytic structure of II n (s) in this model. One can evaluate the Il n (s) on the physical sheet 
by combining Eq.(j5]) and Eq.(j8]) and analytically continue it to the complex plane and to 
the other sheets. Only the poles on the closest Riemann sheet attached to the physical 
region significantly influence the experimental data in a certain physical region between 
two sequential thresholds. For example, in the region between s t h,i and s t h, 2 , one should 
investigate the singularities of the amplitude near this region on the second Riemann sheet. 
In the region between s t h ,2 and s t h, 3 , one should study that of the third sheet, and so forth. 

The picture of this model is clear and easy to understand. In the case that there is only 
one channel (an elastic case), when the bubble chains are summed up, a term (1 —AII(s)) will 
appear in the denominator of the amplitude. When the interaction is attractive and strong 
such that A is larger than l/n(s^), a bound state pole below the threshold appears on the 
real axis of the first Riemann sheet. However, as the coupling becomes weaker, the bound 
state pole will move across the threshold to the second sheet and become a virtual state 
pole. While the interaction is repulsive, there is no bound state pole or narrow resonance 
pole. When more coupled channels are included, these poles usually move to the complex 
s-plane and become resonance poles. 

The experimental data we analyze include the mass distributions of the final states in the 
processes of e + e _ —> X(4260) —> J/^7r + 7r”, (iAZ)*) =l: 7r =F , and . The partial-wave 

decay amplitude could be represented by the formula above with the scattering kinematic 
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FIG. 2. Fitted curves (Solid lines) are compared with the mass distribution data of D°D 


0 n*- 


( D*D*) + represented by the it recoil mass spectrum [6], and J/i^tt [4( from the BES 
group. The dashed line in the first graph represents the contribution of a second-sheet pole at 
(3.875 ± 0.016iGeV) 2 , and the dotted line in the first graph represents that of a bound state pole 
at (3.846GeV) 2 up to normalizations, as mentioned in the text. 


EL 


factors removed. There exist four coupled channels, i.e. r , (DD*^, ( y D*D*) ± , and 
AT(4260)7r , from the lowest threshold to the highest one, and we can define them to be 
channel “1”, “2”, “3”, and “4”, respectively. The X(4260)7 t ± channel is always virtual but 
it provides a background contribution in principle. Since all of the four channels could 
transit into each other through S'-wave, the higher partial waves are assumed to be small so 
that we omit their contributions in the calculation. If every two of them do not couple to 
each other by the s-channel resonance exchange, the effective interactions are represented by 
coupling constants, Ay-s, only ten of which are independent for the time reversal symmetry. 
Care must be taken that in the AT(4260) —* J/iIjtt + 7t~ process , due to the s-t symmetry, 
the invariant amplitude has a contribution from t-ch ann el rescattering as 


M(s ,f, u ) oc Ab(s) + A(A — s — u), 


(9) 


where A = m 2 x ^ 26 ^ + m 2 ^ + 2m 2 , while, in X(4260) —> ( DD *) :f 7r ± and X(4260) —> 
(D*D*) t tt ± processes, only s-channel contributions are considered. 

The total 14 fitted parameters include ten coupling constants, the regularization param¬ 
eter ko, and three normalization factors. The data sets of mass distribution we choose are 
those of J/V'vr - , ( DD *)~, ( D*D*) + (4 6|. The fits to other data sets also leads to similar 
physical conclusions, so we do not show the corresponding fit results here. The best fit 
for the data sets we choose presents a perfect reproduction of the experimental information 
with = 60+3 1 3 2 ^ 4 9 3 _ 11 — 1.03, as shown in Figj2l To reduce the time consumption of 
the numerical calculation, we fixed the k 0 parameter in values between 0.3 and 0.4 in span 
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of 0.01, and found the minimum of the y 2 corresponds to ko = 0.36. The fitted coupling 
constants are listed in the following: 


An = -41.5 ± 0.5, A 12 = -28.0 ± 1.3, A, 3 = 2.5 ± 0.4, 

Ai4 = -12.9 ± 1.7, A 22 = -30.2 ± 4.0, A 23 = -8.6 ± 0.8 

A 24 = -28.6 ± 3.7, A 33 = 71.9 ± 1.1, A 34 = -8.9 ± 1.3 

A 44 = 138.7 ± 1.6. (10) 


With the procedure to extract the pole structure described above, one can easily find the 
nearby pole close to the physical region. The second-sheet poles are searched for near the 
region between s t h,i and s t h, 2 , the third-sheet ones near the region between s t h ,2 and s t h, 3 , 
and the fourth-sheet ones near the region between s^ >3 and s t h,A- A second-sheet pole below 
the DD* threshold and a third-sheet pole just on the DD* threshold are found and the 
locations on the complex s-plane are 

Vs 77 = 3.846 ± 0.019iGeV, 

x/s 777 = 3.875 ± 0.016iGeV. (11) 


However, there is no near-threshold pole close to the D*D* threshold. So the best fit of this 
model prefers a solution with near-threshold poles to produce the Z c (3900) signal, which is 
different from the conclusion of Ref. [12]. 

Until now, there is no method to extract the contribution of every single pole in a multi¬ 
channel amplitude, but we propose that the peak structure in the DD* mass distribution 
is mainly contributed by the third-sheet pole by the following arguments. If the lowest 
J/'ijj'K threshold is omitted, the numbers of Riemann sheet will reduce to 8. Then, the third 
sheet becomes the second sheet and the second sheet becomes the first sheet that the s 11 
is supposed to be mimicked by a real bound state at about (3.846GeV) 2 [18] and s IH a 


second-sheet resonance. I 11 Ref. 


19|, a rigourous formula of unitary scattering amplitude 


with only a virtual state or a bound state in an elastic scattering is shown as 


±- 


T v irtual I bound ( S ) 


s-sl 


sq—sl 

sr—sq 


lTip(s)^yf 


SQ-SL 

SR —SQ 


( 12 ) 


where Sl = (rrq — m 2 ) 2 , Sr = (m 4 + m 2 ) 2 , with m 4 and m 2 denoting the masses of the 
two scattering particles, and s 0 being the mass squared of the virtual/bound state. The 
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amplitude of a resonance in an elastic scattering is 


T 

-L ? 


resonance 


s = 


S G(Zq) 


M 2 (z 0 ) — s — ip(s)sG(zo ) ’ 


(13) 


where 


M 2 (z 0 ) — Re[zo] + 


Im[zo]Im[^op(z 0 )] 

Re[z 0 p(z 0 )] 


G(z 0 ) — 


Im[2 0 ] 


(14) 


Re^op^o)]’ 

and Zq means the pole position on the second sheet. A second-sheet pole at 3.875±0.016GeV 
will contribute a peak abound 3.885GeV, consistent with the experimental data, while a 
bound state at 3.846GeV only contributes a mild enhancement above the threshold as shown 
in Figj 2 ] 

The negative sign of A 22 in this model means a repulsive interaction of DD*, which implies 
that the state might not be formed by DD*. But the D*D* might be attractive for a large 
positive A 33 . By gradually decreasing the interaction strengthes between different channels, 
i.e. the non-diagonal elements of the A matrix, which means that the different channels are 
decoupled with others, the two poles mentioned above move to the same location at real 
axis on different Riemann sheets, yfs = 3.862GeV, which is just the pole position of a bound 
state of D*D* in the elastic case. The pole trajectory means that these two poles are just 
two “shadow” poles originated from the strong attraction between D*D*. 16] It seems that 
the binding energy of the D*D* system is about 140MeV, larger than the usual expected 
value, but it is not in conflict with the heavy quark chiral perturbation calculation. In fact, 
Ref. 20| has exhibited the possibility of binding the vector meson pair though a long-distance 
7 r-exchange potential in heavy quark limit. Furthermore, numerical calculations in Ref. [21] 
show that the bind energy of D*D* through 7 r-exchange could be large up to SlOl^MeV. 

The large positive A 44 also produce a spurious bound state below the threshold 

which could be an artifact of this model, but it does not affect the physical conclusions for 
the reasons below. Since the X(4260)7 t threshold is always virtual in the decay process, this 
channel only contributes to the background of the experimental data. In fact the background 
could come from a lot of other effects such as t-channel contributions and higher partial waves 
and so on. In the fitting process, A 44 is freely running and this channel could absorb all the 
background effects. As a result, A 44 should not be taken seriously as the coupling of this 








FIG. 3. Trajectories of the second-sheet pole and the third-sheet one as all the coupled channel 
effects are being switched off. 


channel. So, it is not surprising that in this channel there could be a faraway pole which 
serves to simulate all the background effects in the physical data region. 

To conclude, utilizing a model summing up all the two-hadron bubble-chain contributions 
and respecting coupled-channel unitarity, we analysed the three invariant mass distributions 
of J/-0Vr 1 * 1 , ( DD *) ± and (D*D*) ± simultaneously in the e + e _ production processes around 
4.260GeV. The numerical result prefers that the Z c (3900) signal to be contributed by a 
pair of near-threshold poles, of which the third-sheet one contribute dominantly. These 
poles might originate from the strong interaction between D*D* through long-distance 7 r- 
exchange interactions, which is similar to the deuteron state. There is no pole structure 
corresponding to the Z c (4025) signal in D*D* mass distribution. Moreover, the unitarized 
scheme may be generalized to analyze other near-threshold structures. 

Note\ When the paper is completed, we noticed a related work is just released in arXiv, 
in which the author also found a second-sheet pole at 3876 + i5MeV in studying elastic DD* 
scattering 22]. 
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